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Abstract. The interaction of a weak probe laser with an inverted-Y type four- level 
atomic system driven by two additional coherent fields is investigated theoretically. 
Under the influence of the coherent coupling fields, the steady-state linear susceptibility 
of the probe laser shows that the system can have single or double electromagnetically 
induced transparency windows depending on the amplitude and the detuning of the 
coupling lasers. The corresponding index of refraction associated with the group 
velocity of the probe laser can be controlled at both transparency windows by the 
coupling fields. The propagation of the probe field can be switched from superluminal 
near the resonance to subluminal on resonance within the single transparency window 
when two coupling lasers are on resonance. This provides a potential application in 
quantum information processing. We propose an atomic ^'^ Rb system for experimental 
observation. 
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1. Introduction 

In recent years, substantial attention has been paid to the study of coherence effects 
in atomic and molecular systems [Ill2ll3llll[5l[6l[71[8l[9]. The interaction of the 
coherent light with multilevel atomic ensembles results in many striking quantum 
phenomena. Three-level atomic systems, such as Lambda, Vee, and cascade schemes 
are the most widely used level schemes for study [TOl [TTl [121 [13] • Among them, the 
electromagnetically induced transparency(EIT), which is based on the phenomenon 
of coherent population trapping [14J, has attracted considerable attention |T5l [161 
[T71 UHl Uni [201 [211 ES]. Some extraordinary effects associated with EIT have 
been studied theoretically, and observed experimentally, including ultraslow pulse 
propagation of light [Il[23], light storage in an atomic vapor [23l[2l], superluminal light 
propagation [HI [25l [26l [271 EHj , coherent control of the optical information processing 
with BEC [29]. An EIT system that results from a quantum interference effect can 
dramatically reduce the group velocity of a propagating probe laser with greatly reduced 
or even vanishing absorption of the probe laser. The essential physical mechanism 
is that the internal structure of atoms can be modified by the interaction with both 
the coupling(or "control") field and the probe field. Therefore the interaction of the 
probe field with the atoms can be manipulated by the coupling field [30]. A properly 
chosen and prepared atomic system is essential for a successful experimental observation. 
Alkali atoms have served as test species for these effects due to the availability of laser 
wavelength and spectroscopy data. Multilevel rubidium atomic systems provide an 
excellent test ground and a starting point for extending the control dimensions with 
inexpensive available laser frequencies for its atomic energy structures. 

In this paper, we investigate the response of a probe laser in an inverted Y-type 
four-level system driven by two additional coherent fields. This scheme has been used 
in the study of the Autler-Townes effect in a sodium dimer [31], and of two-photon 
fluorescence suppression in an ultracold rubidium atom [32], respectively. Here, we 
study the absorption and dispersion of a weak probe laser using probability amplitude 
and equivalent density matrix methods to obtain the linear susceptibility of the probe 
laser. We show that the system exhibits two electromagnetically induced transparency 
windows for the probe field. The transparency windows can be controlled by the 
amplitudes and frequency detunings of the coupling fields. The index of refraction 
associated with the group velocity of the probe laser can be very different at the two 
transparency windows and can be controlled by the two coupling fields. We propose an 
atomic ^^Rb system for experimental observation of this phenomenon. 

2. Equations of Motion 

We consider an inverted-Y type four-level system interacts with three lasers, LI, L2 
and L3 as shown in figure 1(a). Two ground states, |1) and |2), are coupled by laser 
LI and L2 to a common excited state |3), and the excited state |3) is coupled by laser 



Electromagnetically induced transparency in an inverted Y-type four-level system 



3 



L3 to an upper excited state |4), respectively, laser LI is a weak probe laser, L2 and 
L3 are two coupling (or "control") lasers. The corresponding dresscd-statc(|+), |0), 
and |— ))diagram of two coupling lasers interacting with level |2), |3) and |4) is shown 
in figure 1(b). All the transitions are electric dipole allowed. The Hamiltonian of the 
system is given by 

H^Ho + Hi, (1) 

where 

4 

Ho = J2^S){A (2) 

i=l 

is the atomic Hamiltonian, ?uji is the energy of the isolated atom in state and Hj is 
the dipole interaction Hamiltonian, which is given by 

Hi = J2{t\{-fi-E)\j) = -Y^j,,^E,,, (3) 

where /Xjj is the electric dipole moment for \i) ^ \j) transition, and Eij is the 
corresponding coupling laser field. In the rotating-wave-approximation, the interaction 
Hamiltonian can be written as: 

Hi= - ^{Q,e-"'''\3){l\ + Q2e-"'''\3){2\ 

+ Q3e-*^^*|4)(3|) + /i.c., (4) 

where i/j is the laser frequency and fij = ^ijEij/h is the corresponding Rabi frequency 
which is assumed positive in our calculation. The Hamiltonian of the system in the 
interaction representation can be written as 

Hint = - ^(Jlie-^'^i*|3)(l| + Q2e-'''^*|3)(2| 

+ n^e-'''-''\A){'i\) + h.c., (5) 

where 8i = Vi — c<j3i, 82 = 1^2 — ^32, and ^3 = 2/3 — CJ43 are the frequency detunings of the 
probe laser LI, coupling lasers L2, and L3, respectively. uJij = uji — Uj is the \i) ^ \j) 
resonance transition frequency. We assume a;i = for simplicity and the energy of all 
other states are measured relative to state 1 1) . 

2.1. Probability Amplitude Approach 

The atomic wave-function of the system in the interaction picture at any time t can be 
expanded in terms of bare-state eigenvectors as 

l^mt(t)) = ai(t)|l) + a2(t)|2) + a3(t)|3) + a4(t)|4), (6) 

where aj(t) is the time-dependent probability amplitude of the atomic state \i). The 
Schrodinger equation in the interaction picture reads as 

^^ = 4jJ«I*«W> (7) 
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Figure 1. (Color online) (a) Energy level scheme for an inverted Y-type four-level 
atom, (b) Corresponding dressed-state diagram of laser L2 and L3 interacting with 
|2), |3), and |4). 



By introducing the wave-function '^int, and the interaction Hamiltonian Hint of equation 
(5) into the Schrodinger equation, and after making some rotating transformations, we 
obtain the equations for the evolution of probabihty amphtudes of the wave function as 
follows: 

di{t) = i^asi^t) (8a) 
d2it) = ti6i - 62)a2it) + i^asit) (Sb) 
a^(t) = idia^it) + i^aiit) + i^a2{t) + i^aiit) (8c) 

diit) = i{6i + S3)a4{t) + i^as{t) {8d) 

We assume that the probe laser is weak and the population is initially in level 
Oi = 1. We solve the above equations to the first order in terms of the Rabi frequency 
Qi of the probe laser , and to all orders in fl^ and fl2 of the couphng lasers under the 
steady-state condition. From equation (8b)- (8d) we obtain 

4\t) - ^('^1 - S2)aS\t) = 2^a«(t) (9a) 
d«(t) - ^5,a^\t) = 4 + 4a«(t) + ^^a^\t) (95) 
d^\t) - z(5i + 6s)a^\t) = ^^a«(t) (9c) 
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The steady-state solution of a^^ is given by 



The susceptibihty at the probe frequency is given by 

2iV/xi3a*a3 



X 




--T[H^nt,0] + i4).nc, (13) 



di - ^2 ()i + ^3 y 

where N is the atomic number density. Now we include the effects of damping using a 
phenomenological description. Let 73/2 and 74/2 be the decay rates of the probability 
amplitude of levels |3) and |4), respectively. By inspecting equation (9b) and (9c) we can 
see that the effects of damping can be included by replacing Si by Si + ^73/2 in equation 
(9b) and Si + S3 by S1 + S3 + 1^4/2 in equation (9c). This results in the susceptibility as 

^^NjiM^f _ Z73 , , ^^/^ V ...^ 

^ eoh \ ' 2 ^ S1-S2 Si + S3 + i^4/2j ^ ^ 

2.2. Density Matrix Equations 

We can also model this system by density matrix equations. The master equation of 
motion for the density operator in the interaction representation is given by 

where the second term on the left hand-side represents the damping due to spontaneous 
emission and other irreversible processes. In the rotating-wave-approximation it is 
straightforward to obtain the density matrix equations as follows: 

P44 = ^(P34 - P43) - 74P44 (14) 

P33 = ^VPl3 - P3lj + ^(P43 - P34j + -^[P23 " P32) 

+ 74P44 - 73P33 (15) 

P22 = ^{P32 - P23) + VF32P33 - 72P22 (16) 
PU + P22 + P33 + P44 = 1 (17) 
P31 = [iSl - 73l)P31 + ^(Pll - P33) + ^P41 + ^P21 (18) 

P41 = [i{Si + S3) - 74l]p41 + -^P31 - -Y-P43 (19) 

P12 = [i{S2 - Si) - 712JP12 + ^P32 ^Pl3 (20) 
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^ ^^^1 ^^^3 .i^2. . 

P32 = [102 - lZ2)PZ2 + ^Pl2 + ^^42 + ^(P22 " Pssj (^Ij 
P42 = [«((^3 + (^2) - 742]P42 + "^^32 - ^P43 (22) 
P43 = («()3 - 743)P43 ^P41 + ^(P33 - P44) ^p4:2 (23) 

where ji is the population decay rate of level Wij is the branch decay rate 
from level \i) to and •jij represents the coherence decay rate which is given by 
7u ^ 7i« ~ (7* + 7i)/2 + lij- lij the collision dephasing rate. 

We assume that the population is initially in its ground state level |1), and the 
probe laser is weak so that pn (0) ~ 1 for all times. Again, we solve the equations to the 
first order of the Rabi frequency of the probe laser and to all orders of the coupling lasers 
in the steady-state condition. We are interested in the off-diagonal matrix element psi 
associated with the susceptibihty of the probe field. By setting the derivatives to zero 
and keeping the first order term in Qi in equations (18)-(20) we obtain the steady-state 
solution of to the first order of Qi 

(1) _ ^^ipii(o) f , _ , ni/4 , ni/4 



2 V di - ()2 + ^7l2 5i + 53 + i-f 41 ' 



Then the susceptibility reads 



= - ^73l + T /' . + . ^ /' . ■ 25 

eoh V 61-62 + i'ji2 01 + 03 + ^741 J 

The similarity of equation (25) to equation (12) is clear. If the decay rates of the ground 

states(|l) and |2)) and the coUision dephasing rate are small compared to 73, and can 

be neglected, by setting 7^ = 72 = 7?- = 0, equation (25) is identical to equation (12). 

3. Discussion and Numerical Results 

3.1. Absorption Spectra 

It is well-known that the imaginary part of the susceptibility gives the absorption and 
the real part gives the dispersion of the probe field. The susceptibility X equation 
(12) can be separated into the real(x') and imaginary (%") parts as x = x' + ix" ■ The 
explicit expressions for x and x" a-^e 

' _ A^|/^13p('^l-'^2)[A('^l+'^3)+^74/2] . 

„ N\i,,s\''{S,-52)[A^4/2-B{5, + 5s)] 

^ = -7^ ^TW^ ' (2^^) 



with 



+ {S2-S,)M6, + 6,)-^] 
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Figure 2. (Color online) The imaginary part of the linear susceptibility of the probe 
laser as a function of the probe frequency detuning. There are two transparency 
windows at Si = S2, and Si = —S3, respectively, (a) ^2 = 0.673, ^2 = O.573, ^3 = 0, 
and — 73. (b) When 74 = the probe laser at — —63 is also completely 
transparent. Other parameters are the same as in (a), (c) 62 — 0, ^2 — O.573, 
63 — — 0.573, and ^3 = 273. (d) Solid line is the same as in (c); dashed lines are for 
the cascade scheme, |1) — |3) — |4): il2 = and ^3 = — O.573; dotdashed lines are for 
the A scheme, |1) - |3) - |2): O3 = and S2 = 0. 



In our following numerical calculations we assume the atoms are at ultracold 
temperatures and the Doppler effect can be neglected. An ultracold atomic sample 
can be obtained in a magnetic-optical-trap (MOT), such as an ultracold Rb atom trap. 
For example, 10^ ~ 10^ atoms can be trapped in a " dark-spot" MOT within a spherical 
cloud of a size of 1.0 mm in diameter, and the atomic number density N can be around 
10^° ~ [33| A transition dipole moment of 2.5 x 10~^^Cm is used for 

/ii3 in our calculations which is corresponding to the value of ^"^Rb Di transition |35j . 
Combining these parameters, we obtain the coefficient N\pii^\'^ / [tQh) ~ 1.28 x 10^/s. 
For demonstration we assume in our model that the upper state decays much slower 
than the intermediate excited state |3), such as 74 = O.I73 in our calculation. We will 
discuss this with more details later in this paper. 

By inspecting equation (26b), we see that the absorption spectrum of the probe 
laser can have two electromagnetically induced transparency windows as long as the two 
coupling lasers are neither on resonance simultaneously nor at a Raman detuning(52 = 
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— ^3), that is, when the two couphng lasers have different frequency detunings, the 
absorption of the probe laser displays two minima at 61 = 62 and 61 = —63, respectively, 
as demonstrated in figure 2(a)-(c). The absorption of the probe laser becomes zero (or 
the atomic system becomes completely transparent to the probe laser) as the probe laser 
frequency is detuned at Si = 62- The second minima is at 61 = —63 but the absorption is 
not completely zero at this detuning due to the decay of the upper excited state |4). If 
the upper state |4) does not decay, the absorption will also be zero at 61 = —63 as shown 
in figure 2(b). However, the absorption is significantly reduced at 61 = —63 even for a 
decaying upper excited state provided that the coupling laser L3 is strong with respect 
to the decay rate. If we choose an atomic system with a small decay rate 74 of the 
upper state, such as a metastable state, then the absorption of the probe laser can be 
reduced greatly by a relatively strong coupling laser L3 as illustrated in figure 2(c). Let 
us compare this scheme with the widely used three-level A and cascade systems in the 
EIT study. Equation (26) can be used for a A scheme |1) — |3) — |2) by setting fl^ = 0, 
and for a cascade scheme |1) — |3) — |4) by setting fl2 = 0, respectively. In both cases 
there is only one transparency window as illustrated in figure 2(d) by the dotdashed 
line and dashed line. The combination of both systems brings another dimension of 
control of EIT. As one can see from figure 2(d), an absorption line emerges within 
the transparency window of the A scheme by introducing laser L3, and the absorption 
linewidth can be subnatural. 



3.2. The Group Velocity 

It is clear that the absorption spectrum of the probe laser depends on both the detuning 
and the Rabi frequency of the coupling lasers. We plot the imaginary part (absorption) 
and the real part (dispersion) of the linear susceptibility of the probe laser in the same 
frame as a function of the frequency detuning of the probe laser in figure 3(a). The 
dispersion curves(see dashed line) display a very high positive slope at the center of the 
transparency windows; therefore the group velocity of the probe laser can be slowed 
down without absorption. The group velocity of the probe field can be calculated by 
Vg = c/ug, where c is the speed of light in vacuum and the group velocity index is given 

by m 

1 ' 1^1 dx 
2^ 2 dui' 

which is evaluated at the carrier frequency of the probe laser. 

The group index Ug based on equation (27) is simultaneously plotted in figure 3(b) 
as a function of the frequency detuning of the probe laser corresponding to the same 
parameters as in figure 3(a). One can see that the group velocity of the probe laser can 
be reduced by as much as a factor of 10^ for the chosen parameters without absorption. 
When the Rabi frequency of laser L3 increases as shown in figure 3(c), the absorption 
corresponding to the transparency window at 61 = —63 decreases and the width of the 
transparency window increases, while the slope of the dispersion curves decreases at 
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Figure 3. (Color online) Left panel: Absorption (solid line) and dispersion (dashed 
line) of the probe laser (LI); Right panel: the corresponding group index. The 
parameters for the calculations arc: (a)-(b) fls — I.O73, = O.573. (c)-(d) 

= L573, ^2 = 0.573. (e)-(f) r!3 = 1.O73, ^2 = LO73. (g)-(h) = 2.O73, 

V,2 ~ LO73. Other parameters are: 62 — 0, S3 = — O.573, 74 — O.I73. 
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Si = —S3, and therefore the group index at Si — —S3 decreases as shown in figure 3(d). 
When the Rabi frequency of the laser L2 increases the width of the EIT window at 
Si = S2 increases and the central absorption peak is pushed toward to the second EIT 
window at Si = —S3 as shown in figure 3(e), while the group index decreases at Si = S2, 
and increases at Si = —S3 as shown in figure 3(f). When we further increase fl3 the 
central absorption component is pushed back toward the EIT window at Si = S2 and the 
second transparency window becomes wider and deeper as shown in figure 3(g). The 
slopes of the dispersion curves as well as the group index at the EIT windows decreases 
accordingly as shown in figure 3(h). Clearly, the absorption at the two EIT windows, and 
therefore the corresponding group index can be very different depending on the coupling 
lasers. Consequently, the EIT can be controlled by the detunings and Rabi frequencies 
of two coupling lasers as well as the group velocities at these two windows. This can 
be very useful for quantum information processing and transfer. One can control the 
propagation of probe signals at two adjacent frequencies with the two coupling fields. 
To uncover the responsible physical parameters for the group index within each EIT 
window, we write the group index equation (27) explicitly. 

(1) For the Si — S2 transparency window: 

n,{Si = S2) = l + 2Kx(^^^^y (28) 

with K — ^^^Q^^ ; which is a function of S2 and of laser L2, but independent of laser 
L3. For a given detuning S2, the group index at the Si — S2 transparency window is 
inversely proportional to the intensity of laser L2 as illustrated in figure 4(a). 

(2) For the Si — —S3 transparency window: 

n.iSi = -S3) = 1 - 

{Pf-{S2 + S3){a+§^f)) 



q;2 + /32 

«74(^)(52 + 53)/3(a^ + /3f) ^^g^ 



with 



a 



(q;2 + /32)2 



4 

74 



f3 = f - S3{S2 + S3)) 

l^^^^l±3±jn^_s3iS3 + 3S2) 

^ ^ -^3(74 + 73/2) -^.^ 

However, at the transparency window Si = —S3, the group index depends on both 
lasers and it is a function of the three parameters S2, ^2 and Q3 for a given ^3. In figure 
4(b)-(d) we show the group index as a function of fl3, S2, or fl2, respectively, when the 
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Figure 4. (Color online) (a) The group index at the center of the EIT window, 
di — $2, as a function of Q2- Other parameters are 62 — 0.873, S3 — 0, fl^ — 1.673 and 
74 — O.I73. The group index at the center of the EIT window = — ^3 — — O.573: (b) 
as a function of ^l^ for ^2 = and = 0.673; (c) as a function of 62 for il^ — 1.673 
and 0,2 = 0.673; (d) as a function of 0.2, for 62 = and O^ — I.O73. 



two other parameters are given at the transparency window 61 = —63. This clearly 
shows that the group index can be manipulated by the parameters of both coupling 
lasers. 

When both coupling lasers are on resonance or in Raman detuning with each 
other(52 = — ^3), the two transparency windows merge into one as shown in figure 5(a) 
and 5(b), respectively. The width of the transparency window increases with the Rabi 
frequencies, fl2 and ^3, as shown in figure 5(c). If we inspect the spectrum carefully 
we notice that there is a small feature at the center of the transparency window. We 
plot it in an expanded scale in figure 5(d). We find that the dispersion changes between 
normal and abnormal within a very narrow frequency region. Consequently, the group 
index changes from negative to positive; in other words, the probe laser can be switched 
from a anomalous dispersion associated with a negative group index(or superluminal) 
to a positive group index(or subluminal) without being absorbed within this window. 
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Figure 5. (Color online) The susceptibility and the group index as a function 
of the probe detuning when two coupling lasers are on resonance or on Raman 
detuning((52 — —^3)- (a) S2 = S3 = 0, fl^ = — O.573. (b) For laser L2 and L3 
are on Raman detuning: 62 ~ —^z ~ 0.873, Q.-}, — Q.2 — O.573. (c) ^2 — — 0, 
r23 — 2.O73, and ^2 — O.573. (d) we plot (c) on an expanded scale to show the details 
within the EIT window(solid lines: imaginary part of x, dashed lines: the real part of 
X; dotdashed lines: the group index). 



Although the transparency is not one hundred percent for the negative group index 
region the absorption is smalL 

Based on the above analysis, the scheme can be realized in a four-level rubidium 
Eh atom. Two hyperfine ground state levels, |1) = |5S'i/2, F" = 1), and |2) = 155*1/2, 
F" = 2) are coupled by laser LI and laser L2 to a common intermediate excited hyperfine 
level, 1 3) = |5Pi/2, F' = 1), respectively. The third laser L3 couples the intermediate 
level to a higher excited 5D3/2 hyperfine level |3) = \5D3/2,F = 2). In order to observe 
the phenomena experimentally, an ultracold ^^Rb atomic ensemble formed in a MOT 
with a typical temperature of 100 nK would satisfy all the conditions corresponding 
to our above calculations. The lifetime of 51)3/2 state is very long compared to that 
of the 5Pi/2 state(74 ~ O.I73) [35]. So all the conditions assumed in our analysis can 
be satisfied in this system. We expect the experimental observations can be readily 
realized. 
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4. Summary 

In conclusion, we have shown that the response of a probe laser in an inverted Y-type 
four-level system driven by two additional coherent fields exhibits double transparency 
windows for the probe laser. The reliability of the calculations is established by the 
agreement in the susceptibility of the probe laser obtained by both wavefunction and 
density matrix methods. The transparency windows can be controlled by the amplitude 
and frequency detuning of the coupling fields. The group index associated with the 
group velocity of the probe laser can be very different at the two transparency windows; 
hence it can be controlled by the coupling fields. The propagation of the probe field can 
be switched from superluminal near the resonance to subluminal on resonance of the 
probe transition within the single transparency window when the two coupling lasers are 
detuned on resonance. This provides a potential application in quantum information 
processing. This scheme may be realized in an ultracold ^^Rb system and can be used 
to investigate both superluminal and slow light. 
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